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(Dated:) 

We calculate the relativistic entrainment matrix Yik at zero temperature for nucleon-hyperon 

mixture composed of neutrons, protons, A- and E~ -hyperons, as well as of electrons and muons. 

This matrix is analogous to the entrainment matrix (also termed mass-density matrix or Andreev- 
", Bashkin matrix) of non-relativistic theory. It is an important ingredient for modelling the pulsations 

■ of massive neutron stars with superfiuid nucleon-hyperon cores. The calculation is done in the frame 

' of the relativistic Landau Fermi-liquid theory generalized to the case of superfiuid mixtures; the 

£\j i matrix Yik is expressed through the Landau parameters of nucleon-hyperon matter. The results are 

illustrated with a particular example of the o-ui-p mean-field model with scalar self-interactions. 

Using this model we calculate the matrix Yik and the Landau parameters. We also analyze stability 
| of the ground state of nucleon-hyperon matter with respect to small perturbations. 

(vq . PACS numbers: 97.60.Jd, 26.60,+c, 47.37.+q, 97.10.Sj 

<n : 

j~y I. INTRODUCTION 

An analysis of electromagnetic [j], 0, 0, 0] and (in the future) gravitational 0, @, 0] radiation from pulsating neutron 
stars can shed light on the properties of superdense matter in their interiors. The most interesting is the question 
Q-i about the composition of massive neutron-star cores (nucleons? hyperons? quarks? exotic matter?) as well as about 
the properties of superfiuid baryon matter (the dependence of baryon critical temperatures on density, the type of 
5— i ■ pairing of various baryon species). 

' To interpret correctly the observational data it is necessary to have realistic theoretical models of pulsating neutron 
stars. For that, one needs to formulate a hydrodynamics which can be used to describe pulsations. Clearly, the 
ordinary relativistic hydrodynamics (see, e.g., 0]), describing a liquid composed of identical particles, is not suitable 
H ■ f° r this purpose. The neutron-star cores are composed of a mixture of various species of particles with baryons 
J> 1 (nucleons and hyperons) that can be in superfiuid state (0, EH EH El, EH)- The hydrodynamics of superfiuid 
mixtures strongly differs from the ordinary one because it allows the superfiuid components to move independently 
^sD of the normal (non-superfluid) liquid component without any dissipation of energy [lj, EH • 

This paper is devoted to a study of the relativistic entrainment matrix, which is an important quantity in hydrody- 
namics of superfiuid mixtures. We mainly focus on the nucleon-hyperon matter in the core of massive neutron stars. 
Notice that, until now only the superfiuid hydrodynamics of nucleon matter, composed of neutrons (n), protons (p), 
and electrons (e) with a possible admixture of muons (/i), has been considered in astrophysical literature. Let us 
discuss the results of previous works in more detail. 

Assume that the neutrons and protons are in superfiuid state. In this case three independent velocities can exist 
in nucleon matter. Two of them are the velocities V sn and V sp of neutron and proton superfiuid components, 
respectively. The other is the velocity Vq P of 'normal' (non-superfluid) neutrons and protons, as well as electrons and 
muons (it is assumed that, due to collisions, it is the same for all 'normal' particles). The physical meaning of the 
phenomenological superfiuid velocities V sn and V sp can be understood on the basis of microphysics (see, e.g., [16| and 
Sec. IIB). It turns out that the velocity V S i is related to a Cooper pair momentum 2Q i of nucleon species i = n,p by 
the equality 

v si = (i) 

where mi is the mass of a free paricle species i. 

The non-relativistic expressions for the mass current density of neutrons J n and protons J p have the form (see, 

e.g., [H,E3,E1) 

Jn ~ \Pn> ~ Pnn Pnp)Vqp ~\~ PnnVsn ~T~ Pnp^sp? (2) 
Jp (Pp Ppp Ppn)^qp H~ PppVsp PpnV sn- (3) 

Here p n and p p are the neutron and proton density, respectively; pik — ph% is the symmetric 2x2 entrainment matrix, 
also termed Andreev-Bashkin or mass-density matrix (i, k = n, p). It follows from Eqs. ^ and §S§ that superfiuid 
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motion of, for example, neutrons, contributes not only to J n but also to J p (and the same for protons). For the first 
time this effect was predicted, as applied to superfluid solutions of 3 He in 4 He, by Andreev and Bashkin [19]. The 
prefactors in front of V qp in Eqs. ((2|) and ((3]) can be interpreted as the densities of 'normal' neutrons and protons, 
respectively. Since at zero temperature (T = 0) all particles are paired, these densities vanish and we have [lg, [I?], UI| 

Pn Pnn ~r Pnpi (4) 

Pp = Ppp + Ppn- (5) 

More strictly these conditions can be obtained from the requirement of Galilean invariance of the equations of super- 
fluid hydrodynamics at T = [U [H[ . 

The matrix pi k was calculated for the case of T = in Refs. [13, HH and for arbitrary temperatures in Ref. [la |. 
In both cases, the authors used the non-relativistic Fermi-liquid theory of Landau. Though neutrons and especially 
protons in the cores of low-mass neutron stars can be (with a reasonable accuracy) considered as non-relativistic, more 
self-consistent (and necessary in the case of massive neutron stars) is the approach, in which nucleons are treated in 
the frame of relativistic theory. Following Refs. 0, [U HH, [23|, the relativistic analogue of Eqs. @ and §S§ can be 
presented in the form 

3i = I rii - Mfc Y ik U + c 2 ^y ife Q k . (6) 

\ k J k 

Here j i is the particle current density (i = n,p); c is the speed of light; u is the spatial component of the four-velocity 
it M , normalized by the condition u^u^ = — c 2 , and describing the motion of normal part of liquid; rii and [ii are, 
respectively, the number density and relativistic chemical potential of particle species i, measured in the frame where 
U** = (c, 0,0,0). Finally, the symmetric matrix Yi k is the relativistic analogue of the entrainment matrix pi k . In the 
non-relativistic limit Eq. ^ is equivalent to Eqs. © and © under conditions that 

u = V qp , p lk = m l m k c 2 Y lk , (7) 

The prefactor in front of u in Eq. © can be interpreted (by analogy with the non-relativistic case) as the number 
density of 'normal' (non-superfluid) particle species i. At zero temperature this number density vanishes. This imposes 
a condition on the matrix Yi k [2l| 

^PkY ik = rii. (8) 
k 

Taking into account this condition, Eq. © can be rewritten (at T = 0) in the form 

j* = c 2 ]T^Q fe . (9) 

The matrix Yi k for matter composed of neutrons and protons was calculated at T = in Ref. [24| (the authors of Ref. 
[23 | used a somewhat different formalism, see, e.g., the review pBj and references therein). The calculation was done 
in the frame of relativistic a-Lu mean-field model. 

This paper is a natural continuation of the research described above. Our aim is to calculate the relativistic 
entrainment matrix Yi k at zero temperature for matter composed not only of nucleons, electrons, and muons, but 
also of hyperons. We consider only two types of hyperons, namely A- and E~-hyperons (to be denoted by A and S, 
respectively). In most of the calculations, presented in the literature, they appear first in the neutron-star matter 
with the increasing density (see, however, [26|). We wish to emphasize, that analytical results, obtained in this paper, 
can be (in principle) applied to any number of superfluid baryon species. 

At this point it is convenient to make a few remarks concerning the hyperon interactions and superfluidity. First, 
recent experiments indicate that the interaction of E _ -hyperons with nucleons is repulsive (see, e.g., Ref. [26]). For 
a sufficiently strong repulsion it is possible that E~ -hyperons may not appear in the neutron stars at all. In this case 
they can be (in some models) 'replaced' by S~-hyperons (see the discussion of Fig. 11 in Ref. |26|). However, in the 
case of not very strong repulsion, £~-hyperons can appear rather close to a density threshold for A-hyperons (see, 
e.g., Refs. I27L l28ll . The repulsion can also shift the critical temperature of E~-hyperons from (5 x 10 10 — 5 x 10 11 ) K 
(see, e.g., [l3[|29[[3(|) to a lower value. 

Second, as suggested by the Nagara event [3l| , AA attraction can be weaker than it was assumed before. This may 
result in substantial decreasing of critical temperature of A-hyperons (from (10 9 — 10 10 ) K to a temperature much 
less than 10 8 K, see Ref. [l3|). However, it is too early to draw a final conclusion (see, e.g., the criticism of this result 
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on p. 23 of Ref. [13j). The real interaction may be stronger than that deduced from the event (and, of course, new 
experimental evidences are necessary). 

Even if we assume the weak attraction of A-hyperons and the repulsion between E~-hyperons and nucleons, it is 
possible that the inclusion of in-medium effects or three-body forces may significantly influence (increase or further 
decrease) the hyperon critical temperatures. Morever, the exotic pairing of A and S~-hyperons may take place [32l |. 
Taking into account the above discussion, it is reasonable to treat hyperon critical temperatures as free parameters. 
Since in this paper we consider the case of zero temperature, below we assume that all baryon species are superfluid. 

The phenomenological equations ((T|) and ©-([H]), which are discussed in this section in the context of nucleon 
matter, remain unchanged for nucleon-hyperon matter. The only difference is that now the indices i and k run over 
i,k = n, p, A, S (see also Ref. 23]). Thus, now is a 4 x 4 matrix. 

The paper is organized as follows. In Sec. II the relativistic Landau Fermi-liquid theory [33[ is generalized to the 
case of superfluid mixtures. In the frame of this theory we calculate the matrix and express it through the Landau 
parameters f{ k of nucleon-hyperon matter. In Sec. Ill the general results of Sec. II are illustrated with a particular 
example of the a-uj-p mean- field model with scalar self-interactions [34| . Namely, we (i) calculate the matrix Yy.; 
(ii) determine all (spin-averaged) Landau parameters, corresponding to this model; (in) analyze the stability of the 
ground state of nucleon-hyperon matter with respect to small perturbations. Section IV contains a summary of our 
results. 



II. THE RELATIVISTIC ENTRAINMENT MATRIX AT ZERO TEMPERATURE FROM THE 

LANDAU FERMI-LIQUID THEORY 

A. Relativistic Landau theory for mixtures of Fermi liquids 

In this section we briefly discuss how to generalize the Landau Fermi-liquid theory to the case of nucleon-hyperon 
mixture composed of neutrons, protons, A-, and E~-hyperons. The original non-relativistic Landau Fermi-liquid 
theory (e.g., [35|, [3(1]) was extended to the case of mixtures of protons and neutrons in the paper by SjobergjpJ (see 
also [l3|)- The relativistic generalization of the Landau Fermi- liquid theory was given by Baym and Chin |33[ (they 
considered a Fermi- liquid composed of identical particles). 

The generalization of the Landau theory to the case of relativistic mixtures composed of more than one component 
can be made in the same way as in Refs. (33l . [37j . Thus, we only briefly describe the main formulae of the theory 
which will be used subsequently. Notice, that the results obtained below in Section II can be applied to a Fermi-liquid 
composed of any number of baryon species (not necessarily four) . In this case the particle indices in equations should 
run over all baryon species. 

Unless is stated otherwise, throughout the rest of the paper we use the system of units in which the Planck constant 
fi., the Boltzmann constant k^, and the speed of light c equal unity, h = &b = c = 1. We also imply that the subscripts 
i and k refer to baryons. 

As it is demonstrated in Ref. [33|, generally, the structure of the relativistic Landau theory of Fermi liquids is 
the same as that of the non-relativistic theory. Only those results of both theories are different which arc obtained 
using Lorentz (Galilean) transformation properties of various quantities (in particular, the energy and momentum). 
For instance, we will show that the relativistic expression for the effective mass of particle species i differs from its 
non-relativistic analogue. 

Let us consider a system in the ground state with the energy Eq at temperature T = 0. The distribution function 
of quasiparticle species i is then a Fermi sphere, 

n l0 (j>)=0(p Fl -p), (10) 

where p is the quasiparticle momentum; 9(x) is the step function: 6(x) = 1, if x > and otherwise. A small 
deviation Srii (p) of the distribution function from n^o (p) changes the system energy by 

E-E =Y,Zio(p)5ni(p) + ~ fMSmWnkW)- ( n ) 

psi pp' ss' ik 

Addition of one more quasiparticle of a species i, with the momentum p, to the system, increases the total energy E 
by the energy £j(p) of the quasiparticle. From Eq. (fTTj) it follows that 



e,(p) = £io(p) + ]T r k (p,p') Sn k (p'). 

p' s'k 



(12) 
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In Eqs. (fTTj) and (p~2|) p and p' are the particle momenta; s and s' are the spin indices; i, k = n, p, A, and E are 
the baryon species indices. Furthermore, £jo(p) is the energy of a quasiparticle of a species i, corresponding to the 
distribution function n^o (p) • It can be expanded into a series near the Fermi surface in powers of the quantity p — p Fi 
and presented in the linear form 

£io(p) ~ Mi +%0-Ph)> ( 13 ) 

where p H is the Fermi momentum of (quasi)particle species i; fii — Sio(j) Fi ) is the relativistic chemical potential or, 
equivalently, the Fermi energy of quasiparticle species i; v-Fi = [d£io(p) / 'dp\ p=PFi is the velocity of quasiparticles on 
the Fermi surface. It can also be expressed as v Fi = p Fi /m*, where m* is the effective mass of quasiparticle species 
i. Finally, the function / (p,p') in Eq. (fT2"]) is the spin-averaged Landau quasiparticle interaction (here and below 
we disregard the spin-dependence of this interaction since it does not affect our results) . In the vicinity of the Fermi 
surface the arguments of the function p k {j>iP') can be approximately put equal to p » p Fi and p' s» p Fk , while the 
function itself can be expanded into Legendre polynomials Pi (cos 9) , 

f ik (p, P ') = j2fi kp i( cos9 )> ( 14 ) 
i 

where 9 is the angle between p and p' and p k are the (symmetric) Landau parameters, /, = ff l . 

As in the non-relativistic case, the effective mass m* in the relativistic theory can be expressed in terms of the 
Landau parameters /| . To find this relation let us consider, following Ref. [33| . two frames K and K, and assume 
that the frame K moves with the velocity V with respect to K. Below in this section the quantities marked with 
an overline will be referred to the frame K, while those without the overline - to the frame K. The total energy of 
nucleon-hyperon mixture E (E) and its momentum P (P) are related by the Lorentz transformation 

E = (E + PV)j, (15) 
P = P e v (e v P) (1 - 7) + EV 1 . (16) 

In Eqs. (HSJ and (HU) 7 = (1 - V 2 )~ 1/2 , e v is the unit vector along V. 

Now imagine that we add a quasiparticle of a species i, of momentum p and energy £i(p), to the system. Then the 
total momentum and energy in the frame K become equal to P +p and E + £i{p), respectively. On the other hand, 
the momentum and energy in the frame K will be P +p and E + £j(p). Consequently, using Eqs. (p~5|) and (p~6|) one 
obtains the transformation rules for the quasiparticle momentum and energy 

e<(p) = [£,(?>) +i>V] 7, (17) 
p = p-e v (e v p) (1 -7) +£j(p)V7. (18) 

We need to know also how the distribution function of quasiparticle species i transforms from one frame to another. 
The answer is given by the standard formula 

m(p) =ni(p). (19) 

Assume now, that V satisfies the inequality, ^«u H . Then we have also V <C 1 and, as follows from Eqs. (fT7|) — (|T5|) . 
keeping linear terms in V, one gets 

Si(p) = £iO) + ^^£i(p)F-py, (20) 

m(p) = m^ + ^^-e^V. (21) 

In the case of non-interacting relativistic particles the sum of two last terms on the right-hand side of Eq. (|2"D|) equals 
zero, hence £,(p) = £j(p)- 

In addition to Eq. ([20]l . there is one more condition relating £j(p) and e»(p). In fact, it follows from Eq. (|12| that 
for any chosen momentum p the quasiparticle energy ei (p) in the frame if will differ from the energy £j (p) in the 
frame AT only to the extent that rij(p) differs from rij(p). In other words, 



£i(p) = £i(p) + E /**(py) - ^(P')] • (22) 



p's'k 
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Substituting into Eq. (|22"|) e,(p) and nj(p) from Eqs. ([20]) and (j2Tj) , respectively, one obtains 



<9p 



p' s'k 



(23) 



For a system in its ground state one has riiip) — n i0 (p) and £j(p) = £io(p) (see Eqs. (fT0|) and (fl3|) ) . At p = p H Eq. 
([23| relates the effective mass m* and the Landau parameters f{ k 



m* 



(24) 



Here the number density of particle species i is given by 



3tt 2 ' 



while the symmetric matrix Gik equals 



Cm 



— r? r? P k 
g ?r 4 ^Fi^Fk J 1 ' 



(25) 



(26) 



For a liquid composed of identical particles, Eq. (|24p transforms into the equation (13) of Ref. [33|. The non-relativistic 
limit of Eq. (|24p can be obtained if one replaces jij by TOj. Applying then this formula to a mixture of two species of 
baryons, one reproduces the result of Sjoberg [33] (see also jl7l. fl8|). 



B. Calculation of the relativistic entrainment matrix 



Let us employ the theory described above to calculate the relativistic entrainment matrix at zero temperature. 
At first glance, this theory seems inappropriate for calculation of superfluid properties of nucleon-hyperon matter 
because it describes the normal Fermi fluid. However, as was demonstrated by Leggett [38], [39j in the context of 
superfluid 3 He, the particle current density j i of particle species i in superfluid nucleon-hyperon matter is given by 
the same equation as in the case of normal (non-superfluid) matter, namely 



E 

pa 



dp 



m(j>). 



(27) 



All that we need to know is how the superfluid motions modify the distribution function rii(p) of quasiparticles. One 
also has to take into account that a change of rii(p) results in a change of the quasiparticle energy Sj (p). Le ggett 13811 
showed that this energy can be calculated from the same formula (fT2|) as for normal matter (see also [la , [TJIS l24j|). 

As was already mentioned in Sec. I, the superfluid current is generated in the system when the Cooper pairs acquire 
a non-zero momentum 2Q t . In this case they are formed by pairing of quasiparticles with momenta (— p + Q i ) and 
(p + Qi) (rather than with strictly opposite momenta — p andp, as it would be in the system without currents). The 
distribution function rii(p) for a system with currents can be especially easily found at zero temperature. In this case 
all quasiparticles are paired and, up to small terms of the order of 0[(A///) 2 + (Q^/a*) 2 ] (where A is some characteristic 
value of an energy gap in the dispersion relation for baryons; fx is the characteristic chemical potential of baryons), 
rii(p) is a Fermi-sphere, shifted by the vector Q i in momentum space (see, e.g., [171 Il8l. |24|). 



m(p) =0(P K - \P~Q 



(28) 

terms when 



Here and below we assume that Qk <SC p Fi - In this case we may restrict ourselves to a linear in Q k 
calculating j i . Using the distribution function (|28|) as well as Eq. (fl2|) for the energy of quasiparticle species i 
which 8rii(p) = 6>(p H - \p - Q { \) - n i0 (p) w - [dn i0 (p)/dp] Q { , one gets from Eq. (J2ZJ) 



J i 



E 

ps 



dp 



dn M (p) 
dp 



Q, 



E 

ps 



d_ 
dp 



E/ l W)^^Q fc 



p' s'k 



dp' 



(29) 



The first term in the right-hand side of Eq. (|29[) equals / = ni/(m*)Qj. Integrating by parts the second term, one 
has 77 = J2k Gik Qki where the matrix Gik is defined by Eq. (|26|) . Thus, one finds for the particle current density jj 



1 1 1: 



(30) 
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Comparison of this result with Eq. ([9]) allows one to determine the expression for relativistic entrainment matrix Yik 
(5ik is the Kronecker symbol) 

Yik — — r $ik + Gik- (31) 

m* 

Using Eq. we verified that the matrix satisfies the condition ©. 

The energy of nucleon-hyperon matter with superfluid currents can be also expressed through the matrix Yik ■ From 
Eq. ([TTD it follows that 

E-E = ~J2YikQiQ k . (32) 

ik 

An analogous formula, valid for an arbitrary temperature (not only for T = 0) was obtained for a mixture of two 
non-relativistic superfluids by Andreev and Bashkin pj| . Notice that, the difference (E — Eq) can be interpreted as 
the energy of superfluid motion. For a stable superfluid ground state, E — E > 0, and hence the quadratic form in 
the right-hand side of Eq. (|32|) should be positively defined. This leads to a set of conditions on the matrix Yik or, 
equivalently, on the Landau parameters fl k . Here we will write only the simplest two of them (see also (T8L fl9j| ) 

Y u >0, Y u Y kk -Y? k >0 (i^k). (33) 

III. THE RELATIVISTIC ENTRAINMENT MATRIX AT ZERO TEMPERATURE FROM THE o-u-p 

MODEL WITH SCALAR SELF-INTERACTIONS 

Let us apply the general results obtained in Sec. II to a specific model describing the interaction of baryons, the 
a-u>-p mean-field model with scalar self-interactions. Our aim will be to calculate in the frame of this model the 
relativistic entrainment matrix Yik as well as the Landau parameters j\ k of nucleon-hyperon mixture. 

We choose the a-cu-p model rather than, for example, more elaborated mean-field model including hidden strangeness 
a* and <j) mesons (which mediate interaction between hyperons), because of two reasons. First of all, it is relatively 
simple and still realistic model to start with. Second, the hidden strangeness mesons were originally proposed to 
simulate strong hyperon-hyperon interaction. However, Nagara event (3l| suggests that AA interaction in ^^He can 
be weak. To explain such a weak interaction, cr-uj-p model is sufficient (see, e.g., [40|). Bearing this in mind and 
taking into account that the hyperon-meson coupling constants are known with large uncertainty, our choice of the 
model seems justifiable. 

To present a quantitative example supporting our simple model, let us refer to a specific model of neutron star 
cores [27j . This model assumes a weak AA interaction and includes the a* and c6 mesons. As seen in Fig. 6 of Ref. 
[27I ] , the contribution of the a* meson to the hyperon effective masses is at most a few percent of that resulting from 
the a meson. Similarly, the <j) meson potential in which hyperons move is at most a few percent of the contribution 
resulting from the to meson. This example suggests that the contribution of the a* meson to the hyperon entrainment 
matrix is small, while that of the <f> meson may be expected to be small. 

A. a-co-p mean-field model with scalar self-interactions: general equations 

The a-uj-p model with scalar self-interactions is described in detail in the monograph by Glendenning [4l[ (see also 
[3~4j). Here we briefly discuss its main equations which will be used below to calculate the relativistic entrainment 
matrix Yik- Let us consider a system of baryons n, p, A, and £ in some uniform state. Interactions among those 
baryons are mediated by three different kinds of meson fields: scalar cr-field, vector w-field and an isospin triplet of 
charged vector p-fields. The mean-field approximation assumes that the er-, to-, and p-fields are replaced by their 
mean expectation values in the chosen state. We denote these values by er, u>^, and p M = (p^, P2, P3), respectively 
(p, is the space-time index). These mean values are to be calculated from the following (averaged) Euler-Lagrange 
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equations [4l| 



m 2 a a = -bm n g an {ganv) 2 - cg an (c/ CT „cr) 3 

nii - g ai <7 



E 9°i = n< (P) , (34) 

"~ V (P ~ 9W W - ff pl / 3 i P 3 ) 2 + (™i - 5a»Cr) 2 

mf. 



^ = ES . ( 35 ) 



= Pa = 0, (36) 
^E^,jf. (37) 



One sees that only the third isospin component /?g of the p-field, which corresponds to the neutral rho meson, has 
non-zero mean value. In Eqs. (|34[) - (|37p the summation is performed over the baryon species i = n, p, A, and S; mi 
is the mass of meson species I = er, u>, or pi ; 2,3; 9u is the coupling constant of meson I and baryon i; I 3i is the isospin 
projection for baryon species i. Furthermore, u,-(p) is (as in Sect. IIA) a distribution function of particle species i; 
b and c are some dimensionless constants describing self-interaction of the scalar cr-field; w and p 3 are the spatial 
components of four- vectors and pg, respectively. The u>- and p 3 -fields are generated by the baryon four-currents 
j.f on the right-hand side of Eqs. |35|) and (|37jl . They are given by 

fi = n* = 5>*(P)> (38) 



E 



^U(p), (39) 

where the number density n, and the particle current density jj are measured in the laboratory frame; Ei(p) is the 
energy of a baryon species i 



E i(j>) = g^i w + 9 P i hi P3 + V (P _ SU w - P 3 ) 2 + ( m » _ 9oio) 2 - (40) 



In Eqs. (|M|) . I|38|). and (|39|) the summation is performed over the momentum states occupied by the particles. If our 
system is not only uniform but also isotropic then (at zero temperature) the distribution function n.;(p) is a Fermi 
sphere centered at p — in the momentum space, so that we have (see Eq. (|T0|) ) 

n.i(p) = n i0 (p). (41) 

Substituting the distribution function (|41D into Eq. (|38D one obtains that the time component jf = rij is given by Eq. 
(f2l>)) . Moreover, in this special case the spatial components of four- vectors ui^, pg, and vanish, oj = p 3 = j i = 
(there is no preferred direction!), while cr-field and the time components are still given by Eqs. ((34)) . (|35|) . and (|37|) with 
rti(p) and taken from Eqs. (|4"Tj) and J5SJ), respectively. The chemical potential p^ of baryon species i is presented 
in the form 

Mi = &n ^° + Spi ^3* P° + \jP%i + { m i ~ flo-iO - ) 2 - (42) 

It is the energy of a particle on the Fermi surface. 



B. The relativistic entrainment matrix from the o-uj-p mean- field model 

A derivation of the matrix Y^. in the frame of the a-uj-p mean-field model with scalar self-interactions is completely 
analogous to the derivation presented in Sec. IIB for the case of relativistic Landau Fermi-liquid theory. In nucleon- 
hyperon matter in which the supcrfluid currents are generated, the distribution function for baryon species i is 
approximately described by Eq. (|2"5|) . 

The superfluid current density j i is given by Eq. (|39[) with the energy Ei (p) calculated from Eq. (|40| with the help 
of Eqs. $M$>, fH]), and (SJ). 

As it was mentioned in Sec. IIB, we restrict ourselves to a linear approximation when calculating j i as a function 
of Q k . In this approximation the scalar cr-field as well as the time components to and p® remain the same (their 
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FIG. 1: Normalized symmetric matrix Yik/Y as a function of rib for the third equation of state of Glendenning [34]. The 
normalization constant Y = 3no//i n (3no) = 2.48 x 10 41 erg -1 cm -3 . Solid lines show the elements of the matrix Yik/Y; 
each curve is marked by the corresponding symbol ik (i,k = n, p, A, E). Vertical dotted lines indicate the thresholds for the 
appearance of (from left to right) A- and S~-hyperons. 



variation ~ Q^Q^,), whereas the spatial components u> and p 3 depend on some linear combinations of the vectors Q k . 
It follows from Eqs. that 



E 

ps 

E 

ps 

E 



dp 



e{Pn-\p-Qi\) 



dEifr + QJ 



dp 



Mo(p) 



dp 



^Jp 2 + (mi - g al cr) 2 + 



p(Qi - 9uiU - QpihiPs) 



y/p 2 + (mi - g ai cr) 2 



n l0 (p) 



This equation should be supplemented by the expressions ([55]) and l|57|) for to and p 3 , respectively 

Qui A 



UJ 



rn 



2 



V ^ hi 



(43) 

(44) 
(45) 



Solving the system of six equations (J331)~(H3) one can find j t and the vectors w and p 3 as functions of Q k . In this 
way the relativistic entrainment matrix Yn~ can be determined at zero temperature. The analytic expression for Yi k 
is given in Appendix. It is easy to verify that the matrix Y ik satisfies the condition ||5J). 

Note that, in the limiting case considered by Comer and Joynt [H| our results for the matrix Yik do not reproduce 
theirs. Their results do not satisfy the condition (jSJ). Let us remind that the authors of Ref. [24] considered asymmetric 
nuclear matter composed of neutrons, protons, and electrons. They assumed that nucleons interact through a- and 
w-fields (the neutral p3-field and self- interactions of tr-field were neglected). The criticism of such assumption can be 
found in Ref. [42] . 

Fig. 1 presents the normalized elements Yik/Y of symmetric matrix Yik, calculated using Eq. (|rjrj]). as functions 
of the baryon number density rif, — n n + n p + n-s + n\ for the third equation of state of Glendenning [34] . The 
constant Y equals Y = 3no//x„(3n ) = 2.48 x 10 41 erg" 1 cm~ 3 , where no = 0.16 fm~ 3 is the normal nuclear density; 
Mn(3^o) = 1.94 x 10 -3 erg is the neutron chemical potential at n& = 3no- Each curve on the figure is plotted for 
some normalized element of the matrix Yn~ and marked with two particle species indices ik. For instance, symbols 
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FIG. 2: The normalized Landau effective masses m* /mi (i = n, p, A, E) versus n;, for the third equation of state of Ref. [34| . 
Vertical dotted lines indicate thresholds for the appearance of (from left to right) A- and E~-hyperons. 

nA on the figure (right panel) mark a curve plotted for the element Y n ^jY (= Y\ n /Y). The chosen equation of state 
predicts first the appearance of A-hyperons at rib = ribA = 0.310 fm~ 3 and then S _ -hyperons at rib = ribs = 0.319 
fin. . One sees that at rib < nbA (no hyperons) all the components of Yik, related to hyperons, become zero. 



C. Calculation of Landau parameters 



The a-uj-p model described above can be reformulated in terms of the relativistic Landau theory of Fermi liquids 
(see Sec. II). For that, it is necessary to calculate the Landau parameters of nucleon-hyperon matter. In case of 
nucleon matter the Landau par ameters were calculated for various relativistic mean-field models in a series of papers 
(see, e.g., (43L l44l EH. I46L I47L l48l| ) The derivation of these parameters for nucleon-hyperon matter is quite similar. The 
main idea of the derivation is to consider a small deviation of the distribution function of baryon species % from n^o (p) 
(see Eq. [TO)) and to analyze how it modifies the energy of baryon species k. Then the result should be compared 
with the corresponding Eq. (|12p for the energy variation in the frame of the Landau theory. In this way one obtains 
the function f ik (p,p') or, equivalently, the parameters ff k . In Refs. [H, [3, Hi, Hi, H3, IHf , dealing with the case of 
nucleon matter, it is shown that only first two Landau parameters are non-zero: f$ k and f[ k . We checked that the 
same is true for nucleon-hyperon matter, j\ k — at I > 2. In view of this observation it is enough to find only the 
parameters and f[ k . 

Strictly speaking, the parameters f[ k have already been calculated in the previous section. Indeed, it follows from 
Eq. dSJ) that 

9tt 4 / , m 



I'i.i'k v 



fl=Z5Z5-[ Y <k-Z$S*)> ( 46 ) 



where is given by Eq. ([rjrj]) and the Landau effective masses m* (not to be confused with the Dirac effective mass!) 
equal 



- y/p^ + irm-g.icr) 2 . (47) 



\dE l ( J >)/dp\ r _ 



I'Fi 



Fig. 2 illustrates the dependence of normalized Landau effective mass m*/nii (i — n, p, A, S) on nb for the third 
equation of state of Glendenning [34| . 
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Now let us calculate the parameters fr] k - For that we slightly vary the Fermi momentum p H by a small quantity 
Aj3 Fi . This will alter n iQ (p) by 



Sriilp) = 0(p m + Ap H - p) - n i0 (p), 
while the variation of the energy of baryon species i (on the Fermi surface) will be (see Eq. I|12p) 



(48) 



(49) 



Here Srik = p Ffc Ap Fk /ir 2 is the variation of the number density of particle species i. On the other hand, if we consider 
the a-u-p model, the variation of the baryon energy on the Fermi surface will be (in the first approximation, see Eq. 

m) 



SEifaJ = 9uii 5u,° + 9pi I 3i S P ° - Sa, 



(50) 



The small terms 5a, 5uj°, and 8p® can be expressed through Srik from Eqs. (f3"4")l , (|3"5")1 , and (|5T|) . respectively 



Jo, 



L(tr) 



E 



8ni 



E9u)k c 
— on k , 
TY1 



k m p 



(51) 
(52) 
(53) 



The function L(a) in Eq. (|5ip is given by 



5 

i(cr) = [m a a + 6m„ g CTn {g<jn<j) 2 



Cg cm {ganCr) 



E 

psi 



gaijmt - g G io) 
^Jp 2 + {mi - gaio-y 



■ THo(p) 



(54) 



Substituting now Eqs. (I51 j) - ([53j) into Eq. ([50]) and comparing the resulting expression with Eq. (|49|) . one finds the 
Landau parameters fr] k 



f ik _ g^iig^k g P ihi g P khk 
■)a - —3 1 Z3. 



1 



iK - g<,io) gak(m k - g a k<y) 



L(a) 



(55) 



It follows from Eqs. (|4"6"|) and ([55]) that the parameters /q and f{ are indeed symmetric in the indices i and k. 

Just as the parameters f[ k must guarantee the positive definiteness of the quadratic form (f3"2")> . the parameters /o fe 
must satisfy a number of conditions. These conditions are related to stability of charged multi-component mixture 
with respect to density fluctuations and were carefully analysed for nucleon matter (see, e.g., [i^, [5(J HH, HH, HH). 
They depend essentially on the matter composition and on the applied perturbation. Here we consider an equilibrated 
matter of massive neutron stars composed not only of nucleons (n and p) and hyperons (A and E) but also of electrons 
(e) and muons (/i). As an example, we analyse the stability of such matter with respect to long- wavelength density 
fluctuations. 

The stability conditions follow from the requirement of minimum of the free energy F = E — £^ . pjfij (at fixed pj] 
j = n, p, A, E, e, p) for the system in thermodynamic equilibrium, at T = 0. Using Eq. (jlip for the variation of 
energy of baryons, it is easy to find a variation SF = 5E — J^j caused by a small change of Srij(p) (see Eq. P8")l 

with j instead of i): 



"-5E(s*.+J8") 

ik 



Sn^n. + l^iSn^ + ^-p^iSn,) 2 . 



2 On 



2 dn L 



(56) 



Here Ni = m*p H /7r 2 is the density of states of particle species i on the Fermi surface; /i; and n; are, respectively, 
the relativistic chemical potential and number density of electrons (I = e) and muons (I = p). To derive Eq. ([56]) we 
presented the variation SEi of the energy E\ of leptons, in the form (/ = e, p) 



0E1 = — — dm + - -5-5- (dn; ) = /zj dn; + - 7^— (on; ) . 
ojii 2 an, 2 an/ 



(57) 
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same as in Figs. 1 and 2. 



As it should be, the expansion of F begins with the terms of the second order in Srij. The requirement of minimum 
of F means that 8 F > 0, that is the quadratic form in the right-hand side of Eq. (|5rj|) must be positively defined. 

In Eq. (f56"| for the variation SF of the free energy, we neglected a positive term related to the Coulomb energy of 
the perturbed matter. However, it must be taken into account if the perturbed matter acquired a non-zero charge, 
which is the case when 5n p — Sn e — Sn^ — Sn-s ^ 0. The contribution of the Coulomb energy to 5F is then ~ q~ 2 (see, 
e -g-; [49L l50l . I52I]), where q is the wave number of plane- wave density fluctuation. Here we are interested only in the 
limit of long wavelengths, for which q — > 0. In this limit, the positive Coulomb energy can be arbitrarily large, so that 
the matter is stable against the long-wavelength density perturbations at any density. To exclude the 'stabilizing' 
contribution of the Coulomb energy we consider only those variations Snj of the number densities which preserve the 
charge neutrality, 

5n p — Sn e — Sri/j, — Sn-^ = 0. (58) 
Expressing Sn e using this equation and substituting it into Eq. (|56D . one finds 

SF = - ^2 Ajm Srij Sn m , (59) 

jm 

where the indices j and m run over all particle species except for electrons. The 5x5 matrix Aj m is given by 

Aj m = ^-^p + ff^j 5jbS m b + qjq m + 5 m 5 m ^- (60) 

Here Sjb and S m b equal 1 if j and m = n, p, A, £ and otherwise; qj and q m are, respectively, the electric charges of 
particle species j and m in units of proton charge (e.g., q e — —1). 

The requirement of positive dcfinitcncss of the quadratic form (I59[) imposes a set of conditions on the matrix 
elements Aj m or, equivalently, on the parameters /q , we write out only the simplest two of them 

A J3 > 0, (61) 
A 33 A mm - (A jm f > m). (62) 



These conditions are very well known in the literature devoted to stability of nucleon matter (see, e.g., [49l l50l. l52l|). 
For a mixture composed of neutral strongly interacting baryons they can be simplified and presented in the form (see, 
e.g., M) 

1 + i^ 1 > 0, (63) 
(l + i^)(i + ^o fefe )-(^o fc ) 2 >0 (64) 
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FIG. 4: The same as in Fig. 3 but for F{ k . 

where the indices i and k refer to baryons and we introduced the dimensionless Landau parameters Ff , 

F? k = y/NiN k fi k . (65) 

Our results are illustrated in Figs. 3 and 4, where the parameters F^ k and F[ k are presented for the third equation 
of state of Glendenning [34| as functions of rib- The Landau parameters for neutrons and protons are plotted on the 
left panel in Figs. 3 and 4 (i, k— n, p). The right panel demonstrates the Landau parameters related to hyperons 
(i = A, E; k = n, p, A, £). 

We checked that the nucleon-hyperon matter is stable down to baryon number density nj, = 0.34no = 0.055 
fm -3 where the instability occurs (there are no hyperons and muons at such rib). Mathematically, the occurence of 
instability means that the inequality (|62p is not satisfied at < 0.34no = 0.055 fm -3 . Thus, the matter is unstable 
with respect to long- wavelength density fluctuations. All other criteria, which are necessary for positive definiteness 
of the quadratic forms (|32| and (|59|) , are obeyed. 

This type of instability is related to the crust-core phase transition and is carefully analyzed in the neutron-star 
literature (see, e.g., [4!| 0, [H, HI) H3|)- Since we study the stability of matter only in the extreme long- wavelength 
limit and under condition of microscopic charge neutrality, our result for the baryon number density of the crust-core 
interface is just the lower bound for the real value. Precise calculations would give a slightly higher value. For 
example, using extended Thomas-Fermi approach, Cheng et al. [54] found the crust-core boundary at (0.058 — 0.073) 
fm~ 3 , depending on the choice of the a-to-p model parameters. 

IV. SUMMARY 

In this paper we calculated the relativistic entrainment matrix at zero temperature for nucleon-hyperon mixture 
(see Eq. (|3ip ). This matrix is a relativistic analogue of the entrainment matrix (also termed the mass-density 
matrix or Andreev-Bashkin matrix) and is related to pik in the non-relativistic limit by Eq. I|7|. The calculation 
is done in the frame of relativistic Landau Fermi- liquid theory [33|, generalized to the case of mixtures. We show 
that, similarly to pik (see, e.g., [13, EH), the matrix can be expressed through the Landau parameters f\ k of 
nucleon-hyperon matter [i, k = n, p, A, S). If the number of baryon species is more than four, then the indices i and 
k in Eq. <|31|) should run over all these species. 

The general results for Yn,, following from the relativistic Landau Fermi- liquid theory, are illustrated with an 
example of the a-to-p mean-field model with scalar self-interactions. Using this model we obtain the analytic expression 
(1551) for the matrix Y ik . Comparison of this expression with Eq. (|3"Tj) allows to determine the Landau parameters fl k 
corresponding to the chosen mean-field model. Furthermore, we calculate the parameters /q and find that all other 
(spin-averaged) Landau parameters are equal zero, /j = at I > 2. 

In addition, we formulate a number of stability criteria for beta-equilibrated nucleon-hyperon matter (the positive 
definiteness of quadratic forms ([32"]) and l|5"9]) ). Employing the third equation of state of Glendenning 34], which is 
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one of the versions of the a-u-p model with scalar self-interactions, we demonstrate that the nucleon-hyperon matter 
of neutron stars is stable down to the crust-core interface. 

Our results can be used to model the pulsations of cold massive neutron stars with superfluid nucleon-hyperon 
cores. The generalization of these results to the case of finite temperatures will be given in a subsequent publication. 



Appendix A 



Using Eqs. (|43p ~ (|45p . one can express the particle current densities j i as functions of momenta Q k , and thus derive 
the coefficients of relativistic entrainment matrix at zero temperature: 



Yjh — — - 



9wi n k / ffwfc 



4 A ml 



-a-22 



9pkhk 



ai2 



9pihi n k ( g P khk 

A m* k \ m 2 



an 



9uik 
m 2 



-«21 



Here to* is given by Eq. (|4T[) while the coefficients an, ai2, 021, a22, and A are given by 

to 2 to* ' 

9uii9pi^3i n i 



an 



0,12 
ft 21 

O22 

.4 



E 



1 



n 2 J 2 

E9 p i 1 3i ^ 



011022 — ai20.2i- 

In formulae (|67 p - (l70|) the summation is assumed over all baryon species. 



(66) 

(67) 
(68) 
(69) 

(70) 
(71) 
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